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Abstract. We present CCD phase-binned light curves at 
CO 490 nm for 97 Cepheid variable stars in the bar of the 
LMC. The photometry was obtained as part of the French 



in 



EROS project and has excellent phase coverage, permit- 



ting accurate decomposition into Fourier components. We 
identify as 'sinusoidal' or s-Cepheids those stars with pe- 
riods less than 5.5 d and small second-harmonic compo- 
nents. These stars comprise ~30% of our sample and most 
O form a sequence ~1 mag brighter than the LMC clas- 
sical Cepheids in the period-luminosity diagram. They 
5^ are also generally bluer and have lower-amplitude light 
curves. We infer that the s-Cepheids are first-overtone 
& pulsators because, when their periods are converted to 
expected fundamental-mode values, they obey a common 
period-luminosity-colour relation with classical Cepheids. 
This also confirms the reality of the colour term in the 
Cepheid period-luminosity-colour relation. Further, the 
blue edge of the classical Cepheid instability strip agrees 
well with the theoretical calculations for the fundamental 
mode made by Chiosi et al. (1993) for the Hertzsprung- 
Russell and period-luminosity diagrams, but we find that 
our observed s-Cepheids are > 0.2 mag brighter and bluer 
than the Chiosi et al. predictions for the first-overtone. We 
identify a number of features in plots of our stars' Fourier- 
component amplitude ratios and phase differences. These 
features have been identified with resonances between dif- 
ferent pulsation modes. In the LMC we find these features 
seem to occur at periods very similar to Galactic ones for 



classical Cepheids, but at different periods for s-Cepheids. 
We discover a double-mode Cepheid in the LMC, for 
which P(first overtone)/P(fundamental) = 0.710 ±0.001, 
very similar to observed ratios for Galactic double-mode 
Cepheids. 

Key words: Surveys - Cepheids - Magellanic Clouds 



1. Introduction 

Fourier decomposition techniques are proving to be pow- 
erful tools for characterizing the Hertzsprung Progression 
of the changing form of the light curve as a function of 
period exhibited by Cepheid variable stars (Hertzsprung 
1926), and understanding its physical nature and cause. 
For single-mode Cepheids, Simon & Lee (1981) adopted a 
decomposition of the observed magnitudes of form 
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X = X + HfUXi cos ( — i (t - t ) + * 4 ) (1) 

and were the first to point out the existence of considerable 
structure in plots of the amplitude ratios 



~ Xb 
Rki - ^ 



k > 1 



(2) 



and phase differences 
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against period, P, for classical Cepheids in the Galaxy. 
[The form of Eq. (3) renders $ki , which is defined mod- 
ulo 2ir, independent of the epoch, to-] Subsequent work 
- Simon (1988), Poretti (1994), and references therein - 
has shown that with well-determined Fourier components 
there is a clear separation in these plots between clas- 
sical Galactic Cepheids with their skew, high-amplitude 
and possibly bumpy light curves, and another class of 
Cepheids which were recognized by Hertzsprung (1926) 
and which Kholopov et al. in the General Catalogue of 
Variable Stars (1985) define as presenting 'light ampli- 
tudes below 0.5 mag in V (0.7 in B) and almost sym- 
metrical light curves. . . ; as a rule, their periods do not 
exceed 7 days. . . ' In the past some astronomers have be- 
gun refering to these 'symmetrical' or 'sinusoidal' objects 
as s-Cepheids, and we shall adopt this nomenclature, al- 
though a variety of terminologies have been employed. 

Around 1958 Payne-Gaposchkin showed that s- 
Cepheids in the Small Magellanic Cloud are brighter 
than classical Cepheids of the same period (Gas- 
coigne, 1960; Payne-Gaposchkin 1961; Payne-Gaposchkin 
& Gaposchkin 1966) and by 1969 this result had been 
extended to the LMC (Payne-Gaposchkin 1973). It was 
almost immediately suspected that the s-Cepheids might 
be overtone pulsators (e.g. Arp 1960). 

A first mapping of the Hertzsprung progression for 
Galactic classical Cepheids in the R21 — P and $21 — P 
planes was made by Simon & Lee (1981). On the basis of 
a few short-period stars these authors, and Gieren (1982), 
speculated that overtone-mode pulsators might be differ- 
ently located in the R21 — P and $21 — P planes. Then 
Antonello & Poretti (1986), Antonello et al. (1990), Man- 
tegazza & Poretti (1992) and Poretti (1994) determined 
Fourier coefficients for significant numbers of both Galac- 
tic classical and s-Cepheids of short period. They conclude 
that s-Cepheids follow different sequences in the i?2i — P 
and $21 — P diagrams. In the $21 — P plane there is a 
rise in $ 2 i as periods approach ~3.2d. At ~3.2d there is 
a sharp drop in $ 2 i, followed by a further rise at longer 
periods. With a little calligraphic licence we may call this 
form a 'Z-shape'. The drop of the Z-shape is mirrored by 
a minimum in i?2i in the R21 — P plane. It was proposed 
that a separation between classical and s-Cepheids could 
be based on the $21 — P diagram and that the difference 
should be interpreted as the consequence of different pul- 
sation modes: fundamental radial mode ('IF') for classical 
Cepheids, and first-overtone radial mode ('1H') for the s- 
Cepheids. It was further proposed that the s-Cepheid fea- 
tures at P ~ 3.2 d should be interpreted as the signature 
of a resonance between the first and fourth overtone. 

On the other hand, Gieren et al. (1990) believe that 
the overtone nature of Galactic s-Cepheids is established 
only for P < 3.2 d because independent evidence of over- 
tone pulsation has been determined only for some of the 
stars on this branch of the $21 — P diagram. Gieren et al. 



claim that there is clear evidence against the existence of 
s-Cepheid overtone pulsations for P > 3.2 d. 

The differences and similarities between classical and 
s-Cepheids and their pulsational nature can be clarified 
by studying stars in the LMC where observed differences 
in magnitude and colour parallel instrinsic ones. This is 
because the reddening is relatively slight, even in the bar, 
where the IRAS maps of Schwering & Israel (1990) indi- 
cate an uneven distribution of dust. It is also important 
to know to what extent the behaviour of Cepheids in the 
Magellanic Clouds is similar to, or different from, that 
of Galactic ones, since the study of Cepheids at different 
metallicities provides new tests of theories of stellar evolu- 
tion and pulsation. Extragalactic Cepheids are addition- 
ally important because they are acting as a cornerstone in 
establishing the cosmic distance scale. 

Andreasen & Petersen (1987) have studied the Fourier 
components of Cepheids in the Large Magellanic Cloud 
(LMC). Their observational material was 128 photo- 
graphic light curves published by Wayman et al. (1984). 
There were only about 37 points in each B light curve, and 
28 in V . As a consequence of these limitations there were 
only 10 s-Cepheids in their sample and the uncertainties 
on their values of Rki and <J?jti were large. Within these 
constraints, the behaviour of LMC and Galactic Cepheids 
was found to be broadly similar. 

In this paper, we improve upon the Andreasen & Pe- 
tersen results using observations from the EROS project. 
This project (Experience de Recherche d'Objets Sombres, 
Aubourg et al. 1993a, b) has been acquiring photometric 
time series for large number of stars in the LMC in or- 
der to search for possible massive compact objects in the 
halo of the Galaxy via the brightening due to gravita- 
tional microlensing that they should cause of more dis- 
tant LMC stars. However, EROS CCD observations also 
provide high-quality light curves for 97 Cepheids in the 
LMC bar. The high precision and excellent phase cover- 
age (typically 900 points per filter) result in high accuracy 
Fourier components. Further, the fraction of s-Cepheids in 
our sample is much larger (~30%), allowing accurate de- 
lineation of the classical and s-Cepheids in the amplitude- 
ratio and other diagrams. 

2. Observations, reductions and calibrations 

In this paper we use the data from the first season of EROS 
CCD observations in 1991-1992. The EROS CCD equip- 
ment has been described fully by Arnaud et al. (1994a, b). 
The observation, reduction, and calibration procedures 
have been presented in Grison et al. (1995). To summa- 
rize: observations of the LMC bar were obtained at ESO 
La Silla using a 0.4-m, f/10 reflecting telescope and a 
2x8 mosaic of 16 CCDs. Observations were made in two 
bandpasses, and Re, which have mean wavelengths 
of 490 and 670 nm respectively. Pairs of Be and Re im- 
ages were obtained approximately every 25 minutes during 
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clear nights between 1991 December and 1992 April. A to- 
tal of about 1000 images was obtained in each bandpass. 

Because large numbers of stars needed to be reduced, 
a rapid procedure involving the fitting of a Point Spread 
Function was adopted for the extraction of intrumental 
magnitudes at the expense of somewhat reduced photo- 
metric precision. Zero points were established for the nat- 
ural EROS magnitudes Be and Re such that a star of 
zero colour (Be — Re) ~ « (Bj — Vj) has its Be magni- 
tude numerically equal to its Johnson Bj magnitude and 
its Re magnitude numerically equal to its Cousins Rc 
magnitude. In this paper we will mostly use the natural 
magnitudes Be and Re but will also sometimes need to 
convert to or from Johnson Vj and (Bj — Vj) using trans- 
formation equations which can be found in Grison et al. 
(1995): 

Vj = Be- 0A7(B e -Re) o = 0.07 mag (4) 
(Bj-Vj) = 0.92(Be-Re) a = 0.10 mag (5) 

The EROS photometry has a number of limitations. 
We will analyze data only from CCDs 0-4 and 7-10 be- 
cause first-season results from the others were affected by 
a variety of instrumental or reduction problems. The zero 
points of the magnitudes may be systematically in error 
by up to 0.2 mag, and there may be systematic differences 
of the colour zero points of ~0.05 mag between different 
CCDs in the mosaic. Further, crowding or other effects 
can affect the photometry of individual stars significantly. 
In this paper, therefore, it is overall trends which should 
be considered significant, while discrepancies apparent for 
individual stars should only be considered suggestive until 
tested with better- quality photometry with a bigger tele- 
scope. 

3. Identification of Cepheids and Fourier decom- 
position 

In Figure 1, dots show the colour-magnitude diagram for 
stars in one CCD of the EROS mosaic. The main sequence 
and red giant branch are apparent. From the map pre- 
sented by Page & Carruthers (1981) and the article on 
which it is based, we estimate the typical reddening in the 
EROS fields is E(Bj - Vj)=0.15 [« E(B E - Re)] with 
a likely maximum variation of ±0.05. The variations in 
colour and magnitude due to differential reddening are 
therefore very similar to the uncertainties in the corre- 
sponding photometric zero points and cause little addi- 
tional smearing in this or other figures. 

Variable stars were identified using Grison's (1994) 
period-searching procedure, which is able to search effi- 
ciently for periodic, non-sinusoidal light curves by use of 
orthogonal combinations of Fourier harmonics. The limits 
of the Cepheid instability strip were delineated via visual 
examination of a sample of light curves. To avoid contam- 
ination by RR Lyrae stars we only considered variable 
stars with Be < 17.6. Then the light curves of all vari- 
able stars within this region were examined to exclude 



eclipsing binaries. The resulting 97 stars form our sam- 
ple of LMC Cepheids. They are plotted in Figure 1 us- 
ing symbols other than dots. It can be seen that we have 
searched to ~ 1 mag fainter than the lower limit of the re- 
gion where most Cepheid or Cepheid-like stars occur. Our 
Cepheids are listed in Table 1. The naming and numbering 
scheme continues that introduced by Grison et al. (1995). 
The J2000 equatorial coordinates were derived via Har- 
vard (x, y) coordinates in the manner described by these 
authors, and are accurate to ~3 arcsec. 

Mean light curves were obtained by sorting the ~900 
best photometric points in each colour into phase order 
and then taking the mean magnitude averages in bins 
of 0.01 in phase. Figure 2 presents the Be light curves, 
which are less noisy than the red ones. Only 29 of our 
Cepheids were previously identified by Payne-Gaposchkin 
(1971) on Harvard plates (HV stars). Comparison of our 
periods with those quoted by Payne-Gaposchkin suggests 
our typical period uncertainty is about 0.005 days. 

Coefficients derived from the Grison period-searching 
procedure were used to calculate the parameters Xi, $j, 
Rki and for various values of M, the order of the 
Fourier decomposition (defined in Eq. 1). Tests showed 
that the parameters were only negligibly affected by 
the value of M (for M greater than 5), except for the 
few bump Cepheids in our sample (EROS 2081, 2087, 
2089...). 

We adopted a procedure whereby M was increased 
until there ceased to be any decrease in the variance a 
of the residuals between the mean light curves and their 
Fourier models. Typical values of a are in the range 10~ 2 
to 5 10~ 3 mag. Values of R21, R31, $21 and $31 for the 
Be light curves are listed in Table 1. Their uncertainties 
were calculated as described in the Appendix. Also listed 
in Table 1 are magnitude means < Be > and < Re > and 
the peak-to-peak amplitudes ABe and A_Re- 

We can give an estimate of the Cepheid detection ef- 
ficiency, or, equivalently, the lowest amplitude detectable. 
These can be derived from the work of Groth (1975) 
and Scargle (1982). From the statistical properties of the 
periodogram if a signal is present in the data (Groth 
1975), Scargle (1982) has defined a detection efficiency 
which gives the probability of detecting a periodic signal, 
acos(wi), for a given signal to noise ratio, a/a, where a 
is the amplitude of the signal, and a the average preci- 
sion of the measurements. Applying such results to the 
EROS data we find with a confidence level better than 
99% that we expect to detect all the periodic sinusoidal 
variable stars with a/a ~ 0.3. The average precision a 
of the ~900 CCD measurements has been evaluated from 
the least-squares residuals of the Cepheids' light curves. 
Thus, the theoretical detection limit of a sinusoid with 
a peak-to-peak amplitude A = 2a is typically 0.01 mag 
for a Cepheid of 14th mag and 0.02 mag for one at 16th 
mag These predictions are in agreement with the lowest- 
amplitude contact eclipsing binary found in the EROS 
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Fig. 1. The colour-magnitude diagram of mean values for all stars reduced from CCD 1 for which < Be >< 18.5 (dots). The 
entire sample of Cepheids is also shown: Classical light curves (diagonal crosses), s-Cepheids (filled diamonds), double-mode 
Cepheid (open circle), intermediate objects in the R21 — P plane (asterisks) and anomalous types (open triangles). The same 
symbols are used in subsequent figures. The full, dotted and dashed lines represent the blue edges of the instability strip for 
fundamental, first and second overtone pulsation, transformed from the calculations by Chiosi et al. (1993). 



data (EROS 1003, A ~ 0.04, B E = 14.7; Grison et al. 
1995). We have not estimated our detection efficiency for 
double- mode Cepheids. 

4. Classification 

4-1. Anomalous stars and double-mode candidates. 

The distinction between the skew, higher-amplitude light 
curves of classical Cepheids and the more-sinusoidal, lower 
amplitude ones of s-Cepheids can be clearly seen in Fig- 
ure 2. Pairs of (Classical Cepheid, s-Cepheid) with nearly 
identical periods include (2019, 2020), (2022, 2023), (2034, 
2035) and (2066, 2067). However some light curves appear 
unusual and are likely to confuse the discussion of the sep- 
aration between classical and s-Cepheids. 

EROS 2037 and 2076 have markedly larger scatter in 
their mean light curves than other stars of similar period. 
This might be due to photometric contamination by su- 
perposed or nearby stars. Another possibility is that the 
stars might be double- mode Cepheids. 

We searched for a second periodicity in these light 
curves. We subtracted the Grison Fourier model obtained 
on the mean light curve from the unbinned complete data 
and then reapplied the Grison algorithm to the residu- 
als. Only EROS 2037 showed a second period, of value 
2.443 ± O.OOld. For this star the ratio of the second to 
first periods is then 0.710 ±0.001, in close agreement with 
the observations reviewed by Balona (1984) for Galac- 
tic double-mode Cepheids. This strongly suggests that 



this star, for which irregular variations were first noted 
by Shapley & McKibben Nail (1955), is a double-mode 
Cepheid. It has also been classified by Alcock et al. (1995) 
as a double-mode Cepheid pulsating in the fundamental- 
mode and first overtone. The value of the period ratio they 
obtain is 0.710, like the value we derived. 

Andreasen & Petersen (1987) list 7 candidate LMC 
double-mode Cepheids, but none has been confirmed. Al- 
cock et al. (1995), using data from the MACHO microlens- 
ing project, report the detection of 45 double-mode stars 
amongst 1500 Cepheids. We find only one secure double- 
mode variable out of 72 stars with P < 5d. 

The Fourier coefficients for EROS 2026, 2027, 2051 and 
2075 may be affected by significant systematic uncertain- 
ties because of the gap in the light curves and possible 
period errors as a consequence of these stars' periods be- 
ing close to an integral number of days. Some examples 
of unsatisfactory phase coverage and its effect on Fourier 
decomposition can be found in Antonello et al. (1990). 

Additionally, two of the stars already identified as 
unusual (2095 and 2097) fall clearly below the general 
trend in the magnitude-period diagram plotted in Fig- 
ure 4. Probably they are W Virginis stars associated with 
the LMC according to the calibration of Demers & Har- 
ris (1974). In all the figures in this paper, unusual objects 
are plotted as open triangles, with the exception of the 
double-mode star 2037, which is plotted as an open circle. 



J. P. Beaulieu et al.: EROS Cepheids in the bar of the LMC 



LU 
Q 
ZD 



CD 
< 



LU 
DC 
LU 



1 1 1 1 A 1 1 1 1 1 1 1 

: 2001 "\. x : 
i.307d : 


: 2002 ♦ : 

. 1 .364d ■ 


: 2003 ♦ : 

1.531d '- 


: 2004 ♦ : 
1 .620d : 


: 2005 ♦ : 

- / \ - 

\ 

1 .664d '- 


: 2006 ♦ : 

W \*.w ***** *"\v 

1 .736d - 


: 2007 ♦ : 
. .- 

1.765d 


I I I I I I I I I I I 
: 2008 ^ ♦ : 

V...-' 

: 1 .905d 


I I I I I I I I I I I 
: 2009 ♦ : 

1 .906d 


I I I I I 
: 2010 ♦ 

i.9i4d : 


I I I I I I I I I I 
2011 ♦ : 

2.080d : 


I I I I I I I I I I I 
: 2012 ♦ : 

X 

2.193d : 


: 2013 ♦ : 

2.242d 


I I I I I I I I I I I 
: 2014 ♦ : 

: ..•-•""' v x.. : 

2.344d 

I I I I I I I 


I I I I I I I I I I I 
: 2015 ♦ : 

*-"\-'..V 

2.41 2d '- 

I I I I I I I I I I I 


I I I I I 
: 2016 ♦ : 

2.41 8d : 

I I I I I 


I I I I I I I I I I 
: 2017 ♦ : 

." r "N. 

\ 

t / \ - 
2.481 d '- 

I I I I I I I I I I 


I I I I I I I I I I 
: 2018 ♦ : 

:""-.~-'" v 

2.558d : 

I I I I I I I I I I 


: 2019 .*\ x : 

2.620d 


: 2020 ♦ : 

2.630d ' 

I I I I I I I I I I I 


: 2021 x : 

■-. 

"■"•.v.- -■' * v -v^ 
-.V 

2.730d - 

I I I I I I I I I I I 


: 2022 ,. x x : 

2.737d ' 


: 2023 ♦ : 

2.866d " 

I I I I I I I I I I 


: 2024 ♦ : 

2.951 d '- 

I I I I I I I I I I 


: 2025 ' ••, a : 

' "'"--/ 2.963d 


: 2026 a : 

2.976d ' 

I I I I I I I I I I I 


: 2027 \ a : 

2.993d " 

I I I I I I I I I I I 


: 2028 .' \ x ; 

\ 

3.088d - 


: 2029 ♦ : 

3.088d " 

I I I I I I I I I I 


: 2030 x : 

..--\ 

3.103d '- 

I I I I I I I I I I 


: 2031 ;' \ x : 

3.246d : 


: 2032 ,-. x : 

\ 

"*.». * 

3.266d : 


: 2033 ..\ x : 

3.281d '- 


ri 

: 2034 .-' x : 

3.367d : 


: 2035 ♦ : 

3.394d '- 


: 2036 ♦ : 

,. / \ 

3.399d : 


: 2037 •. o : 

3.444d : 


I I I I I I I I I I I 

; 2038 x : 
""""^ 

: 3.482d : 


I I I I I I I I I I I 

: 2039 ♦ : 

3.519d '- 


I I I I I I I I I I 
'. 2040 * . 

3.535d : 


I I I I I I I I I I 

: 2041 a x : 

3.550d ! 


I I I I I I I I I I 
: 2042 ; \ x : 

3.567d - 


: 2043 ♦ : 

3.598d : 


I I I I I I I I I I I 

2044 _„ x : 

\ 

3.604d : 


I I I I I A I I I I I 

: 2045 •' \ x : 

"W N 

3.657d - 


'. 2046 ■•' x : 

3.71 9d - 


I I I I I I I I I I 

: 2047 /\ x : 

3.795d - 


I I I I I I I I I I 

'. 2048 x : 

- 

3.820d : 


: 2049 ♦ 
3.863d : 




I I I I I I I I I I I 

: 2050 X 

\ 

"\ 

3.874d : 


I I I I I I I I I I I 

2051 \ a : 

3.983d '- 


1 1 1 I II I 1 I 1 I 
2052 / \ x : 

4.051 d : 


I I I I I I I I I I 
2053 x : 

s. 

4.1 08d : 


I I I I I \ I I I I 

2054 ;'\ X : 
• 

4.113d: 



-0.3 


+0.3 

-0.3 


+0.3 

-0.3 


+0.3 

-0.3 


+0.3 

-0.3 


+0.3 

-0.3 


+0.3 



0.5 



0.5 



0.5 



0.5 



0.5 



PHASE 



Fig. 2. EROS differential Be light curves for 97 Cepheids in the bar of the LMC averaged in 0.01 bins of phase. The value of 
the magnitude mean < Be > is reported in Table 1. The curves have been phased such that the rising branch attains < Be > 
at phase zero. 
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PHASE 

Fig. 3. (a) The fundamental mode for EROS 2037 
(P = 3.44 d). The points show the full, unbinned Be 
light curve. The full line shows the Fourier model to the 
phase-binned data of Fig. 2. 



4-2. Classical, sinusoidal and indeterminate Cepheids 

Figures 5 and 6 show plots of i?2i and $21 against period. 
It is evident that the points fall in several groups, and as 
noted above, earlier authors have adopted various classi- 
fication criteria to distinguish classical from s-Cepheids. 

It seems to us that the most transparent discrimination 
on the form of the light curve should be in the i?2i — P 
plane where a low value of -R21 indicates more sinusoidal 
light curves. Mindful of earlier work on Galactic Cepheids, 
we have selected as s-Cepheids those stars lying in the 
region defined by i? 2 i < 0.30 for P < 3d and i? 2 i < 
0.20 for 3.0 < P < 5.5 d. [This last upper limit is not 
well-defined by our data. Work on Galactic Cepheids (e.g. 
Simon & Lee 1981) shows that the i?2i value for classical 
Cepheids - as defined by the period-luminosity relation - 
falls to low values for P ss 7 — 10 d.] Our s-Cepheids are 
plotted in all figures as filled diamonds. Our sample is 
composed of 30 s-Cepheids according to this criterion. 
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Fig. 4. Period-luminosity relation for the 97 EROS Cepheids in the Be filter (mean wavelength=490 nm; same symbols as Fig. 
1). Most s-Cepheids are ~1 mag brighter than classical Cepheids of the same period. The full, dotted and dashed lines represent 
the blue edges of the instability strip for fundamental, first and second overtone pulsation, as calculated by Chiosi et al. (1993) 
with no convective overshoot. 




5. Discussion 



5.1. Pulsation modes 



Fig. 3. (b) The overtone mode for EROS 2037 plotted at the 
same vertical scale as Fig. 3a. The points represent the mean 
Be light curve in 0.01 bins of overtone-mode phase after sub- 
traction of the Fourier model obtained for the fundamental 
mode. The zero of phase is arbitrary. The full line represents 
the Fourier model for the overtone mode taken to the third 
order. The ratio of periods P(1H)/P(1F) = 0.710 ± 0.001. 



A group of 6 stars with J? 2 i « 0.24-0.27 and 4 < P < 
6.5 d (2040, 2057, 2061, 2065, 2072 and 2078) stands apart 
as of uncertain appartenance in Figure 5. They are plotted 
as asterisks. 

We classify the remaining 51 stars as classical Cepheids 
and plot them as diagonal crosses. 



For the vast majority of stars, we see in Figure 4 that 
our morphological classification based on Fourier param- 
eters is mirrored by a clear separation into two period- 
luminosity sequences, with the s-Cepheids ~1.0 mag 
brighter than the classical ones at any given period. Figure 
6 shows the classification is further mirrored by a separa- 
tion in the $21 — P diagram, and we shall see below that 
this is true of other Fourier parameters too. 

An obvious interpretation is that the lower sequence in 
Figure 4 corresponds to fundamental-mode pulsators and 
most classical Cepheids (as we define them) are pulsat- 
ing in this mode, while the upper sequence corresponds to 
first-overtone pulsators and most of our s-Cepheids pul- 
sate in this mode. Figure 1 then shows that the overtone 
pulsators are in general bluer, as suspected for Galactic 
Cepheids by Antonello et al. (1992). 

The few stars which deviate from this correspondence 
clearly merit further observation to determine whether 
they really do deviate, or whether their EROS photometry 
is systematically in error, as discussed in Section 2. Exam- 
ples are EROS 2023, and possibly EROS 2010, which have 
s-Cepheid values of R21 but fundamental-mode values of 
luminosity and colour, and EROS 2033, 2055 and 2070 
where classical Cepheids appear to have first-overtone lu- 
minosity and colour. Examination of these latter three 
stars' images shows that they are oval. Thus they have al- 
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Fig. 5. The R21 — P diagram for Be light curves for the 94 objects with P < 14 d (same symbols as Fig. 1). We define as 
s-Cepheids those non-anomalous stars within the region defined by R21 < 0.30 for P < 3d and R21 < 0.20 for 3.0 < P < 5.5 d. 
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Fig. 6. The $21 — P diagram for Be light curves for the 94 objects with P < 14 d. The classical and s-Cepheids follow different 
loci. 



most certainly been brightened by optical binarity, which 
could also explain their bluer colours. 

Figure 4 indicates that the intermediate or indetermi- 
nate stars with _R 2 i ~ 0.24-0.27 are a mixed bag: five 
are fundamental-mode pulsators, one (2057) is an over- 
tone pulsator, and one (2037) is pulsating in both modes 
simultaneously. Figure 7 confirms that in general the s- 
Cepheids have lower amplitudes than the classical ones 



and shows that the intermediate stars, besides having in- 
termediate values of R21, have intermediate amplitudes. 

The identification of s-Cepheids as overtone pulsators 
is supported by the theoretical models of Chiosi et al. 
(1993). In Figures 1 and 4 we have drawn in the blue edges 
of the fundamental and first and second overtone insta- 
bility strips transformed from models with no convective 
overshoot (which Chiosi et al. claim best-matches LMC 
Cepheids), Y=0.25, Z=0.008 and the Padova scale. In 
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Fig. 7. Amplitude-period diagram for the 94 objects with P < 14 d. The s-Cepheids have a strong tendency to lower amplitudes 
than the classical Cepheids. 



both figures the fundamental-mode blue edge delimits the 
distribution of classical Cepheids well. The double-mode 
Cepheid EROS 2037 even falls close to the fundamental- 
mode blue edge. However the EROS observations show 
that first-harmonic blue edge is both brighter and bluer 
than the Chiosi et al. predictions by > 0.2 mag. 

More-detailed comparisons are hardly justified given 
the zero-point uncertainties in the EROS photometry, but 
we note that Chiosi et al. models with mild, and especially 
full overshoot (not plotted) bound the observed distribu- 
tion of classical Cepheids in Figure 4 less well. 

5.2. 'Modified period'- Luminosity relation. 

We have plotted the Figure 4 period-luminosity relation in 
a different way in Figure 8. We have modified the period 
of the presumed overtone pulsators to the corresponding 
fundamental period, using the ratio P(1H)/P(1F) « 0.70. 
The classical and s-Cepheid sequences are now closer, 
but not coincident. This is because a fundamental and 
first-overtone Cepheid of the same magnitude are not 
otherwise-identical stars differing only in period: they also 
differ in colour, as seen in Figure 1. 

5.3. 'Modified period' -Luminosity-Colour relation. 

Sandage (1958) showed that Cepheids are more accurately 
represented by a period-luminosity-colour relation. Nu- 
merous attempts have been made to determine the re- 
lation. It is not critical which one we use: we have chosen 
the period-luminosity-colour relation for LMC Cepheids 



determined by Feast (1984): 

<M V >= alog(P) +(3{<B > - <V >) + (6) 

with a=-3.8, /3=2.70 and </>=-2.39. 

In Figure 9 we take account of the colour term. The 
period plotted for a given star was determined in the same 
manner as in Figure 8. The ordinate is < Vj > — f3(< Bj > 
— < Vj >) + const, obtained using Eqs. 4 and 5. If the 
PLC relation is applicable, the expected slope in Figure 9 
is a. 

The diagram has considerably lower scatter than Fig- 
ures 4 and 8, and we remark that most of the classi- 
cal and s-Cepheids lie on the same sequence. The colour 
term in the PLC relation is clearly justified, and the same 
colour correction applies to first- overtone and to funda- 
mental pulsators. 

Some stars do not follow the relation. This may be due 
to optical or physical binarity. Another possibility is that 
the star is a second-overtone pulsator (e.g. 2001, 2003, 
2015, 2067). 

5.4- Fourier parameters 

5.4.1. The $21 — P plane: fundamental pulsators 

The classical Cepheids, or fundamental pulsators, fall 
along a strip corresponding to the Hertzprung progression. 
From comparison with Galactic results (Simon 1988, and 
references therein ) we would expect $21 to rise sharply 
to 2ir at P ss 10 d and thereafter rise from zero towards a 
constant value for longer periods. This steep change has 
been interpreted as a 2:1 resonance between the second 
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Fig. 8. Period-luminosity relation. The periods used in the diagram for the first-overtone candidates and first-overtone blue 
edge were modified using the relation P/0.70. Compare with Fig. 4. 
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Fig. 9. Period-luminosity-colour relation. The periods were modified as in Figure 8. The ordinate is the numerical value 
< Be > — 2.95(< Be > — < Re >) + 2.39. Compare with Figures 4 and 8. The classical and s-Cepheids are seen to obey the 
same PLC relation for the fundamental period. 



overtone and the fundamental in the normal mode spec- 
trum of these stars (Simon 1988, and references therein). 
Andreasen & Petersen (1987) and Andreasen (1988) have 
shown that in both the LMC and the SMC the resonance 
occurs around 10 days. From our sample we can conclude 
that the resonance occurs between 8.7 d (lower limit given 
by EROS 2091) and 10.2 d (upper limit given by EROS 



2092). A better value will be given by a bigger sample of 
LMC Cepheids. 

For LMC Cepheids with periods greater than 12 days, 
EROS 2093, 2094, and 2096 suggest that the limiting value 
of $21 is about 4rad, which is similar to the values found 
by Andreasen & Petersen (1987) for 7 stars with periods 
greater than 12 days observed photographically. 
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5.4.2. The $21 — P plane: overtone pulsators 

For overtone pulsators, two well-defined sequences noted 
by Antonello et al. (1990) appear. The upper sequence is 
characterized by P < 3.2 d, $21 > 4.2 rad, and for the 
lower sequence 3.0 < P < 6d and $21 < 4.0 rad. These 
upper and the lower sequences, and their possible link is 
what we called the Z-shape. It is of great importance to 
notice that the stars on the horizontals of the Z near where 
they join the diagonal are located at two or more times 
their formal error from the classical sequence. Therefore 
the existence of the Z-shape in the $21 —P plane is beyond 
doubt. 

From our earlier discussion we can see that the stars 
belonging to the lower part of the Z, the lower sequence, 
have high luminosity and are overtone pulsators. There- 
fore, the alternative hypothesis of Gieren et al. (1990) that 
these stars could be fundamental pulsators is refuted. 

The drop in the Z-shape is observed in the Galaxy at 
a period around 3.2 d and has been interpreted as the sig- 
nature of a possible 2:1 resonance between the first and 
fourth overtone. A comparison of the Z-shapes exhibited 
in the Galaxy and the LMC gives information about the 
effect of metallicity on the position of this possible reso- 
nance. 

The position of the drop is difficult to estimate. In our 
sample of s-Cepheids, we decided not to take account of 
EROS 2023, which despite its low R21 would appear to be 
a IF pulsator, because its magnitude is faint, and EROS 
2024, because of the large uncertainties on its $21 value. 

Visual inspection of the diagram places the drop in 
the range 2.5-3.4 days. In the $21 — P plane, the limit of 
the lower sequence is given by EROS 2029 (P = 3.0882 d, 
$21 =3.38 rad), and a limit on the upper sequence can be 
found from EROS 2018 (P = 2.5582 d, $21 =4.85 rad). We 
do not observe a period overlap between the upper and the 
lower sequence. EROS 2020 (P=2.6303d, $ 2 i =3.83 rad) 
falls between these limits, on top of the classical sequence, 
about 3 standard deviations below the upper sequence. In 
our opinion, its observed $21 value is a signature of the 
joining branch between the upper and the lower parts of 
the Z-shape. With only one object in this branch it is not 
possible to give a certain value for the position of the drop. 
We propose that it may be at P « 2.7 ± 0.2 d, that is, at 
a shorter period than in the Galaxy. Buchler & Moskalik 
(1995) observed that in the SMC the drop should occur 
at a lower period than in the Milky Way, but were unable 
to give a precise position for it because they only have a 
few objects in the lower part of their Z-shape. 

The EROS CCD data from the 1992-1994 observa- 
tional campaigns should double the number of Cepheids 
and lead to a better definition of the overtone sequence. 



5.4.3. The R 31 - P, $31 - P and $41 - P planes. 

In Figures 10, 11 and 12 we present R31 — P, $31 - P 
and $41 — P diagrams. Classical-Cepheid features in these 
planes have already been discussed by earlier authors. Si- 
mon & Moffet (1985) have shown that some resonances 
may present stronger signatures using $31 and $41 rather 
than only $ 2 i- A sharp dip feature appears in _R 2 i — P 
at P ss 7.5 d, though its appearance is confused by the 
slower fall of R21 towards the minimum at P ~ 10 d. The 
dip is possibly weakly visible in $21- By going to upper 
harmonics in the $31 — P and $41 — P planes, the dip ap- 
pears more and more clearly. This sharp dip involves five 
EROS stars (2083-2088). It is centered at 7.5 ± 0.2 days. 

Moskalik et al. (1992) used the Iglesias-Rogers opaci- 
ties in order to study the beat and bump Cepheid mass 
discrepancy via non- linear calculations. They found a dip 
in the Rzi — P plane at 7.7 d and related it to the sharp 
feature observed (but not interpreted) by Simon & Moffet 
(1983) in the $41 — P plane near 7 days. Moskalik et al. 
consider that it can possibly be related to the near 3:1 res- 
onance between the fundamental and the fourth overtone. 
We conclude that our dip is a possible signature of this 
3:1 resonance observed for the first time in the LMC, at 
the same period as the feature observed in the models of 
Moskalik et al. 



6. Conclusion 

We have Fourier analyzed the light curves of Cepheids in 
the bar of the LMC using photometry obtained by the 
EROS project. This is the first study of a large sample of 
Cepheids in the LMC using high-quality photometry ac- 
quired over a long time scale with a large CCD mosaic. 
Good phase coverage was obtained for almost all the ob- 
jects. 

Using these light curves we studied Cepheid pulsations 
in a nearby galaxy with a metallicity different from our 
own. Because of the low differential reddening and known 
distance of the LMC, differences in apparent magnitude 
correspond to differences in absolute magnitude, and it 
was possible to derive new and convincing results con- 
cerning fundamental and overtone pulsations previously 
only suggested by comparable studies in the Milky Way 
by Poretti (1994), Simon (1988) and others. 

The colour-magnitude diagram was used to select 97 
Cepheid or cepheid- like variable stars. Only 29 of them 
were previously recognised in Payne-Gaposchkin's (1971) 
Harvard Catalogue. 

Eight of the stars present unusual light curves, and 
one is a double-mode Cepheid. Almost all the rest divide 
cleanly into two groups in the R21 — P diagram. We as- 
sociate the group with periods less than about 5.5 d and 
low values of i?2i with the s-Cepheids and the rest with 
classical Cepheids, though this definition is not of course 
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Fig. 10. The R31 — P diagram for Be light curves for the 94 objects with P < 14 d. The sharp dip for classical Cepheids at 
P « 7.5 d is stronger than in Fig. 6. 
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Fig. 11. The $31 — P diagram for Be light curves for the 94 objects with P < 14 d. The sharp dip for classical Cepheids at 
P « 7.5 d is stronger than in Fig. 6. 



phenomenologically identical to that adopted by other au- 
thors. 

The stars in the two groups mostly follow separate, 
parallel loci in the period-luminosity diagram, with the 
s-Cepheids being more luminous by ~lmag for any given 
period. It has been suggested that s-Cepheids are first- 
overtone pulsators, and theoretically it is expected that 
the ratio of periods P(1H)/P(1F) ~ 0.7. When the s- 
Cepheid periods are modified by this factor and a stan- 



dard LMC period-colour-luminosity relation is applied, 
the classical and s-Cepheid loci are found to be conin- 
cident, furnishing persuasive evidence that s-Cepheids are 
pulsating in the first-overtone and that the colour term is 
justified in the PLC relation. Further, the overtone pul- 
sators are found in the colour-magnitude diagram to oc- 
cupy the bluer portion of the instability strip, indicating 
that it is colour which determines whether a Cepheid pul- 
sates in its fundamental or first harmonic. The blue edge 
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Fig. 12. The $41 — P diagram for Be light curves for the 94 objects with P < 14 days. The sharp dip for classical Cepheids 
at P « 7.5 d is stronger than in Figures 6 and 11. 



of the fundamental-mode instability strip agrees well with 
the theoretical calculations of Chiosi et al. (1993), but 
the first-overtone blue edge is observed to be bluer and 
brighter by > 0.2 mag than the Chiosi et al. predictions. 

The overtone pulsators also tend to have smaller am- 
plitudes and are separated from fundamental pulsators in 
the $21 — P diagram. 

Seven stars have intermediate values of R21 and light- 
curve amplitude. They prove to be of varied nature: 5 
fundamental-mode pulsators, one overtone pulsator and 
one double-mode Cepheid. 

We identify about 50 fundamental and 30 obvious 
overtone pulsators, and most of these stars have periods 
less than 10 d. Payne-Gaposchkin & Gaposchkin (1971) 
showed that the LMC was overabundant in short-period 
Cepheids compared to the Galaxy, but detected only 3% 
of s-Cepheids. This low proportion undoubtedly reflects 
the difficulty of detecting low amplitudes in photographic 
data. We find, as suspected by Matteo et al. (1990) on the 
basis of a sample of 6 stars, that the proportion of LMC 
overtone pulsators is much higher, representing ~30% of 
our sample, and ~40% of stars with P < 5.5 d. 

Some Cepheids are brighter than expected from the 
modified period-colour-luminosity relation. It is known 
that the EROS photometry can be uncertain for individual 
stars, but the deviant objects might also be Cepheids with 
optical or physical binary companions, or possibly cases 
of second overtone pulsators. These possibilities could be 
tested with better-quality photometry. 

Other stars requiring further photometry include 
EROS 2023, and possibly EROS 2010, which have sinu- 
soidal light curves but classical-Cepheid luminosities. If 



these luminosities are confirmed, it will indicate that IF 
pulsators can occasionally have symmetric, low-amplitude 
light curves. 

The features exhibited in the R21 — P and $21 — P 
planes are similar to those found for Galactic Cepheids. 
Most of our stars have period less than 10 days: however, 
it is possible to constrain the position of the 2:1 resonance 
between the fundamental and second overtone for classical 
Cepheids. We think that this resonance occurs between 8.7 
and 10 d, as already suggested by Andreasen & Petersen 
(1987). 

A sharp feature near 7.5 d is observed in the Fourier 
planes (especially the higher-order ones). We suggest it 
is related to a similar feature observed in the models of 
Moskalik et al. (1992), which in turn is believed related to 
the near-3:l resonance between fundamental and fourth 
overtone. 

The Z-shape around ~3 days in the $21 — P diagram 
is clearly seen. The overtone nature of s-Cepheids with 
period greater than 3 days is confirmed by their higher 
luminosity than classical Cepheids. Moreover, the 'drop' of 
the Z-shape for LMC first-overtone pulsators is estimated 
to be at ss 2.7 ± 0.2 d. This value is smaller than the one 
observed for Galactic Cepheids. This may be due to the 
lower metallicity of the LMC. 

We thus find that the position of resonances in over- 
tone pulsators is shifted in the LMC, whereas for funda- 
mental pulsators the values are similar to Galactic ones. 

Using the 1993-1994 EROS CCD data, we expect to 
double the number of Cepheids in our sample. This will 
better define the positions of the resonances. Moreover, 
comparison with observations of similar nature carried out 



18 



in the SMC, like the work of Buchler & Moskalik (1994), 
should be of great interest. 
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Appendix. Uncertainties for the Fourier decompo- 
sition parameters 

For a model of the form given by Eq. (1), uncertainties for 
the Fourier decomposition parameters R k i and have 
already been calculated and discussed by Petersen (1986). 
In this paper we take a different model, because the Gen- 
eralized Periodogram method introduced by Grison (1994) 
requires a different form of Fourier decomposition. 

Grison's procedure is to take a model for a curve of 
zero mean and of form 
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ized Periodogram method. The aj and bj are random inde- 
pendent variables. Thus, they are uncorrelated. Therefore, 
by approximating the function Rki = f(cii, bj) by its first- 
order Taylor development around < Rki >, the variance 
o\ can be given by the expression: 



where 



R k i 



9 Rki 

da,i 

dRki 
dh 



M 

E 



9 Rki 

dcii 



ol + 



9 Rki 



(8) 



9R k i 



and 



Rki 



da* 



9Nr, 



( D Rkl \ 1 1 


— Rki 


1 


\N Ru , ) D Ru , 


-ft-fcl 


(N Rhl y i 


1 


1 


\D Rk J D Rkl 


Rki 


DR k i 


( 




i < k 


< 2A k 




i = k 



dD 



R k i 



dai 

9D Rk 
9h 



-2( a ° ki A k + 0° ki B k ) 


2B k 

-2(a' ki A k + (3 s ki B k ) 
2Ai 

-2{aliAi+^ u Bi) 
2Bi 



> k 
< k 



> k 
= 1 

> 1 
= 1 

i > 1. 



-2(aliAi + ^Bi) 

The same reasoning applies to the phase differences $ k i 
for which it is easy to show that 



$i = arctan ( 



-Bi 



V Ai 



1,...,M. 



(9) 



(We define over [— tt, tt] ; the sign of is given by the 
sign of —Bi.) Therefore the variance o% is given by the 
relation: 



M 



E 



where 



da 

9§ki d$k 



°l + 



9$ki 



(10) 



with 



dN h 



with 



2 ( f Rki 



N, 



9cii 
9N Rkl 



9h 



9R k i 



9R k i 



dD, 



dai 
dD Rkl 



dh 



--{Al+BD 
D Rhl = (Al + B\) 



9$k 

dai 



dai 



9$k 
dh 



dai 
dh 


h k B k 



9at 
d$ k 
9h 



9<&i 
dai 
gji 
dh 



i < k 
i = k 



^(-alfik+^Ak) i>k 



J^i-a^Bi + pc u Ai) i>l 



i < k 
i = k 

(-a s ki B k + f3 s ki A k ) i>k 



h k 
A k 



J. P. Beaulieu et al.: EROS Cepheids in the bar of the LMC 



' l^(-"f<Si+/9f^i) i>l. 

and 

ht = 5 • 

References 

Alcock C, Allsman R.A., Axelrod T.S. et al., 1995, AJ in press 
Andreasen G.K., 1988, A&A 196, 159 
Andreasen G.K., Petersen J.O., 1987, A&A 180, 129 
Antonello E., Aikawa T., 1993, A&A 279, 119 
Antonello E., Poretti E., 1986, A&A 169, 149 
Antonello, E., Poretti, E., Reduzzi, L., 1990, A&A 236, 138 
Antonello E., Poretti E., Reduzzi L., 1992, Mem. Soc. Astron. 
Ital. 63, 75 

Arnaud, M. Aubourg E., Bareyre P., et al., 1994a, Experimen- 
tal Astron. 4, 265 

Arnaud, M. Aubourg E., Bareyre P., et al., 1994b, Experimen- 
tal Astron. 4, 279 

Arp, H.C., 1960, AJ 65, 404 

Aubourg E., Bareyre P., Brehin S., et al., 1993a, The Messenger 
72, 20 

Aubourg E., Bareyre P., Brehin S., et al., 1993b, Nat 365, 623 
Balona L.A., 1984, The double-mode Cepheids. In: Madore, 
B.F. (ed) Proc. IAU Colloq. 82, Cepheids, theory and ob- 
servations. Cambridge University Press, Cambridge p. 17 
Buchler J.R., Moskalik P., 1994, A&A 292, 450 
Chiosi, C, Wood, PR., Capitano, N., 1993, ApJS 86, 541 
Feast, M.W., 1984, Magellanic Cloud Cepheids. In: van den 
Bergh, S., de Boer, K.S. (eds) Proc. IAU Symp. 108, Struc- 
ture and evolution of the Magellanic Clouds. Reidel, Dor- 
drecht p. 157 
Gascoigne, S.C.B., 1960, Trans. IAU 10, 680 
Gieren W.P., 1982, PASP 94, 960 

Gieren W.P., Moffet T.J., Barnes T.G., Matthews, J.M., Frueh, 

M.L., 1990, AJ 99, 1196 
Grison P., 1994, A&A 289, 404 

Grison P., Beaulieu J. P., Pritchard J.D., et al., 1995, A&AS 
109, 1 

Groth, E.J., 1975, ApJS 29, 285 

Hertzsprung, E., 1926, Bull. Astron. Inst. Netherlands 3, 115 

Iben, I, Jr., Tuggle, R.S., 1975, ApJ 197, 39 

Kholopov P.N., Samus' N.N., Frolov M.S., et al., 1985, Gen- 
eral Catalogue of Variable Stars, 4th edition, Vol. I. Nauka, 
Moscow 

Mantegazza I., Poretti E., 1992, A&A 261, 137 
Mateo M., Olszewski E.W., Madore B.F., 1990, ApJ 353, Lll 
Moskalik P., Buchler J.R., Marom M., 1992, ApJ 385, 685 
Page, T., Carruthers, G.R., 1981, ApJ 248, 906 
Payne-Gaposchkin, C.H., 1961, Vistas Astron. 4, 184 
Payne-Gaposchkin, C.H., 1971, Smithson. Contr. Astrophys. 
No. 13 

Payne-Gaposchkin, C.H., 1973, Comparison of the Cepheid 
Variables in the Magellanic Clouds and the Galaxy. In: 
Muller, A.B. (ed) The Magellanic Clouds. Reidel, Dordrecht 
p.34 

Payne-Gaposchkin, C.H., Gaposchkin, S., 1966, Vistas Astron. 
*8, 191 



19 

Petersen J.O., 1986, A&A 170, 59 
Poretti E., 1994, A&A 285, 524 
Sandage A., 1958, ApJ 127, 513 
Scargle, J.D., 1982, ApJ 263, 835 

Schwering, P.B.W., Isreal, F.P., 1990, Atlas and catalogue of 
infrared sources in the Magellanic Clouds. Kluwer, Dor- 
drecht 

Shapley, H., McKibben Nail, V., 1955, Proc. Nat. Acad. Sci. 
41, 185 

Simon N.R., 1985, ApJ 299, 723 

Simon N.R., 1988, The Fourier decomposition technique: Inter- 
preting the variations of pulsating stars. In: Stalio, R., Will- 
son, L.A. (eds) Pulsation and Mass Loss in Stars. Kluwer, 
Dordrecht p. 27 

Simon N.R., Lee A., 1981, ApJ 248, 291 

Simon N.R., Moffett, T.J., 1983, PASP 97, 1078 

Stobie, R.S., 1977, MNRAS 180, 631 

Wayman PA., Stift M.J., Buchler C.J., 1984, A&AS 56, 169 
Westerlund B.E., 1993, The stellar populations in the Mag- 
ellanic Clouds, an overview. In: Baschek, B., Klare, G., 
Lequeux, J. (eds) New aspects of Magellanic Cloud Re- 
search. Springer- Verlag, Berlin p. 7 



This article was processed by the author using Springer- Verlag 
MEN A&A style file L-AA version 3. 



